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Abstract 



Two compounds, tetragonal YB4 and cubic YBg, have been investigated by electric- 
field gradient (EFG) and Knight shift measurements at the boron sites using the ^^B 
nuclear magnetic resonance (NMR) technique and by performing first-principles cal- 
culations. In YBg ^^B (/ = 3/2) NMR spectra reveal patterns typical for an axially 
symmetric field gradient with a quadrupole coupling frequency of uq = 600ibl5 kHz. 
In the second boride (YB4) three different EFGs were observed corresponding to the 
three inequivalent crystallographic sites for the boron atoms {Ah, 4e, and 8j). They 
correspond to: j^q(4/i) = 700ib30 kHz with an asymmetry parameter r] = 0.02ib0.02, 
z/Q(4e) = 515 ± 30 kHz, t] = 0.00 + 0.02/ - 0.00, and ugiSj) = 515 ± 40 kHz, 
r) = 0.46 ± 0.08. The Knight shifts measured by Magic-Angle Spinning (MAS) 
NMR at room temperature are very small being 0.6 it 8 ppm and —1 it 8 ppm for 
YB4 and YBe, respectively. For the theoretical calculations structure optimizations 
were performed as a first step. For the obtained structural parameters the EFGs 
were computed within the local-density approximation. Very satisfactory agreement 
between experimental and theoretical results is obtained both for the structural pa- 
rameters and the B EFGs thus confirming the underlying structural models. In 
addition to the EFGs, band structures, densities of states, and valence-electron 
densities are presented and the bonding situation in the two yttrium borides is dis- 
cussed. The band-structure results are compatible with the very low values for the 
Knight shifts mentioned above. 
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1 Introduction 



Through their specific properties, transition metal borides became attractive 
both from scientific and technical point of view. For instance, YBe, like some 
other rare-earth and actinide hexaborides, shows superconductivity below 
7.1 K [1] and has been considered as a possible candidate for high-temperature 
thermoelectric materials [2]. 

Detailed knowledge of the electronic structure of a material is a key ingredient 
for an in-depth understanding of many of its macroscopic features. A combi- 
nation of experimental methods as well as computations on an ab-initio level 
has proved very efficient, in particular relating electronic structure calcula- 
tions to results of nuclear magnetic resonance (NMR) technique. The value of 
the electric field gradient (EFG), measured by NMR quadrupole interaction, 
is directly determined by the charge distribution around the nucleus. Thus, 
theoretical EFG studies based on the electronic structure are important in 
order to give a reliable interpretation of the experimental data. One of the 
strengths of NMR measurements is that, being a microscopic tool, they arc 
sensitive to the symmetry of the crystallographic sites and in particular to the 
electron density in the vicinity of the nucleus. In a recent study, this sensitivity 
was used to provide criteria for the determination of crystal structures that 
are still under debate [3,4,5]. 

In the present paper, we report the EFG values for YBg and YB4 and the 
determination of the Knight shifts by using the Magic- Angle Spinning (MAS) 
technique. The experimental EFGs are compared and interpreted with the 
theoretical ones obtained from ab-initio calculations. 

The structure of cubic YBg (space group PmSm, no. 221) consists of a simple 
cubic lattice of corner-connected Bg octahedra with Y atoms fiUing the cubic 
holes. In this and other borides with the same structure the B-B distance be- 
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tween neighbouring octahedra is slightly shorter than within the octahedron 
(for YBg: 1.64 A compared to 1.75 A). In agreement with this crystallographic 
structure there is only a single EFG value. A previous ^■'^B NMR study for 
YBq [6] reported the spin-lattice relaxation time and the spectrum. However, 
the authors were unable to determine the EFG because of the widely broad- 
ened sateUite lines and the weak NMR signal used in their instrumentation. 

On the other hand, tetragonal YB4 (space group PA/mhm, no. 127) has a 
quasilayer structure with alternating sheets of Y atoms and Bg octahedra 
linked together laterally by B2 units. In this structure five different nearest- 
neighbour B-B distances between 1.64 A and 1.81 A are observed. In YB4 
the boron atoms are situated in three different crystallographic positions with 
different EFGs for each of them, as will be demonstrated in the present paper. 



2 Structure optimization 



For confirmation of the validity of the available experimental lattice param- 
eters and the atomic positions, calculations have been performed using the 
Vienna ab-initio simulation package (VASP) [7,8,9]. By this method the Kohn- 
Sham equations of density-functional theory [10,11] with periodic boundary 
conditions are solved within a plane-wave basis set with electron-ion inter- 
actions described by the projector augmented wave (PAW) method [12,13]. 
For exchange and correlation the general gradient approximation (GGA) [14] 
was applied. The structural parameters were calculated by atomic forces and 
stress-tensor minimization. For each of the two borides two different energy 
cutoffs for the plane-wave basis were used (400 cV and a cutoff higher by at 
least a factor of 2), which lead to practically the same results, so that adequate 
convergence is ensured. The experimental and calculated structural parame- 
ters for YB4 and YBg are given in Tables 1 and 2, respectively. For cubic YBg, 
where the lattice parameter a and the positional parameter of the B atom are 
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Table 1 

Experimental and calculated structural parameters for YB4 adopting the ThB4 
structure {P4/mbm, no. 127). The lattice parameters are in A. 



Lattice 


Positional 








parameters 


parameters 


Atom 


Site 


Remarks 


a = 7.111 


X = 0.3179 


Y 


45 


exp., Ref. [15] 


c = 4.017 


z = 0.2027 


B(l) 


4e 






X = 0.0871 


B(2) 


4h 






X = 0.1757 


B(3) 


8j 






y = 0.0389 








a = 7.1035 








exp., Ref. [16] 


c = 4.0206 










a = 7.1091 


X = 0.3182 


Y 


45 


calc, VASP (GGA) 


c = 4.0280 


z = 0.2030 


B(l) 


4e 






X = 0.0870 


B(2) 


4h 






X = 0.1760 


B(3) 


8j 






y = 0.0386 









the only free parameters, the minimum-energy structure was also found by 
the full-potential linearized augmented plane- wave (FLAPW) method within 
the LDA approximation which is described in the next Section. 

The structural parameters optimized by VASP have been used to calculate 
the EFGs for YBg and YB4. 
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Table 2 

Experimental and calculated structural parameters for YBg adopting the CaBg 
structure {Pm3m, no. 221). The lattice parameters are in A, the positional param- 



ctcr refers to site (3/ of B. 


Lattice 


Positional 




parameter a 


parameter x 


Remarks 


4.102 




exp., Ref. [17] 


4.1025 


0.202 


exp., Ref. [18] 


4.113 




exp., Ref. [19] 




0.199 


exp., Ref. [20] 


4.0436 


0.1989 


calc, FLAPW (LDA) 


4.1004 


0.1988 


calc, VASP (GGA) 



3 Electric-field gradients 

3. 1 Experimental 

To avoid skin-depth effects for better RF penetration, the samples were used 
in powder form and obtained by crushing the single crystals. The crystals 
were obtained by crucible free inductive zone melting of the YB4 and the YBg 
source rods that provides a high purity and single-phase material. They, in 
turn, were prepared in the following steps: 

(1) Mixing of Y2O3 (purity 99.999%) and boron components in proper ratios. 

(2) Synthesis of the YB4 and the YBe powders by routine solid-state reaction 
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of borothermic reduction: 



Y2O3 + IIB ^ 2YB4 + 3B0 t 

Y2O3 + 15B ^ 2YB6 + 3B0 T 

(3) Slip casting of source rods, removal of the binder and sintering of the 
rods in the vacuum. 

The ^^B NMR experiments were carried out with a Bruker DSX Avance spec- 
trometer at a frequency of 96.29 MHz. The static spectra were obtained by the 
Fourier transform of the free induction decay (FID) following a short single 
pulse (in the range of 0.9-2.5 /is). The spectra contain up to 256 accumulations 
with a repetition time of 20 s. A 4 mm MAS probe was used with rotation 
frequencies between 7 and 9 kHz. The chemical shifts are given with respect 
to external BFgEtaO. In Fig. 1 the ^^B NMR spectra for YB4 and YBg are 
given. 

3.2 First-principles calculations 

The all-electron band-structure calculations for the calculation of EFGs are 
based on the density-functional theory [10,11] (DFT) and the local-density 
approximation and have been performed by the linearized augmented plane- 
wave (LAPW) method [21] in its full-potential version [22,23,24,25] (FLAPW) 
using an exchange-correlation potential by Hedin and Lundqvist [26,27]. 

The EFGs have been calculated from the I — 2 components of the Coulomb 
potential near the nuclei. The formalism by Herzig [28] and Blaha et al. [29] has 
also been employed to split the calculated EFG components into the contribu- 
tions from the surrounding electrons within the respective muffin-tin sphere 
( "sphere contribution" ) and the remainder that comes from outside this sphere 
("lattice contribution"). This partitioning depends, to a small extent, on the 
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choice of the muffin-tin radii. The valence contribution can be spht further 
into the allowed W contributions (only sd, pp, pf are important in the present 
context) which provide useful information about the influence of particular 
Z-like wave functions on the EFGs [30]. As is common practice the EFG com- 
ponent with the largest absolute value is always designated as Vzz- For further 
algorithmic details on the EFG calculation see our recent paper [4]. 

3.3 EFGs: Results and discussion 

Figure 1 shows the ^^B spectra which are typical for a nuclear spin / = 3/2 in 
the presence of first-order quadrupole effects. For YBg (bottom of Fig. 1) the 
separation of the satellite lines is given by Su — i/q (3 cos^ — 1), where © are 
the powder singularities. © = 90° thus yields uq — e'^qQ/2h. Here, eg = Vzz 
and Q is the largest component of the electric-field-gradient tensor and the 
nuclear quadrupole moment, respectively. The axially symmetric field gradient 
is consistent with the local symmetry at the boron site. The existence of three 
inequivalent boron sites in YB4 is the origin of the more complex satellite spec- 
trum compared to YBg (Fig. 1). The i/q value and the asymmetry parameter 
V — iVxx~Vyy) l^zz wcrc obtained using the dmfit simulation program [31] . The 
vq value translates into Cq = 2uq and it follows for Vzz = 4.136 x lO^^C^/Q 
(in V/m^) when the appropriate Q value (in barns) for ^^B is used and Cq is 
inserted in units of MHz. Since different nuclear quadrupole moments can be 
found in the literature, the experimental values for Vzz were determined for 
two different values, namely Q = 0.04 b and Q — 0.0355 b [32]. The former 
Q value has been taken from the almanac of the spectrometer manufacturer 
Bruker and is very close to Q = 0.04059 [33] recommended by Pyykko [34] in 
his compilation of nuclear quadrupole moments. The results are shown in Ta- 
bles 3 and 4 for YBg and YB4, respectively. Experimentally only the absolute 
value of Vzz and 77 can be determined for the powder specimens. 
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Table 3 

Calculated B EFGs (in 10^° V/m^) for YBg for the optimized structure compared 
to the experimental results assuming an ^^B nuclear quadrupole moment of 0.04 
and 0.0355 |e| x 10~^^ m^, respectively. 



Calc. 


Exp. 


Exp. 


Opt. struct. 


Q = 0.04 


Q = 0.0355 








-13.5 


12.4 ±0.3 


14.0 ± 0.3 



Table 4 

Calculated B EFGs (in lO^'^ V/m^) for YB4 for the experimental structure [15] 
and the optimized structure compared to the experimental results assuming an ^^B 
nuclear quadrupole moment of 0.04 and 0.0355 |e| x 10~^^ m^, respectively. 







Calc. 


Calc. 


Exp. 


Exp. 






Ex]:). struct. 


Op1. s( 


ruct. 


Q = 0.01 


Q = 0.0355 


Site 


i 






Vii 


V 


l^nl V 




B(l) 


z 


-11.0 


0.0 


-10.8 


0.0 


10.6 ±0.6 0.0 


12.0 ±0.6 0.0 


B(2) 


X 


-8.1 




-8.2 










y 


-8.3 




-8.4 










z 


16.4 


0.02 


16.6 


0.01 


14.5 ± 0.6 0.02 


16.3 ± 0.6 0.02 


B(3) 


X 


2.7 




2.9 










y 


7.9 




7.9 












-lO.G 


0.48 


-10.8 


0.17 


10.6 ± 0.9 0.46 ± 0.09 


12.0 ±0.9 0.16 ±0.08 



For YBe the B EFG corresponds to the axially symmetric case (site symmetry 
C^y) where there is only one independent EFG component. In YB4 the situa- 
tion is more comphcated. The three crystallographically distinct B sites (site 
symmetries: C4 for B(l), C2V for B(2), and Cg for B(3)) lead to one, two, and 
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Table 5 

Split of the B EFG for YBq (optimized structure) into lattice and sphere components 
and the latter into its main contributions, i. e., sd, pp, and pf. All Vzz values are in 
units of 10^° V/m^. 



Site 


^zz ^ ZZ 


^ zz 


^ zz 


T/PP T/P/ 
^ zz " zz 


B 


-13.5 9.0 


-22.5 


-1.9 


-17.8 -2.3 



three independent EFG components, respectively. This means that for B(2) 
and B(3) a non-zero asymmetry parameter 77 is observed. 

The comparison between calculated and measured EFGs shows very good 
agreement, if the absolute values are considered. In general for the EFG tensor 
either one or two components are negative. The corresponding principal axes 
are determined by the strongest bonding interactions. In the case of YB4 the 
negative EFG component for B(l) {Yzz) is related to the axis of the inter- 
octahedral B-B bond along the c axis, for B(2) the two negative components 
iYxx and Vyy) refer to the (001) plane through B(2) and its three B neighbours, 
and finally for B(3) the negative EFG component (Vzz) refers to the axis of 
the bond between the adjacent atoms of a Bg octahedron and a B2 unit. 

Now the splits of the B EFGs into the lattice and sphere contributions and the 
latter into their main components {sd, pp, and pf) are considered (Tables 5 
and 6). The sphere contribution is the one with the largest absolute value, 
but also the lattice contribution is relatively large. This behaviour has been 
observed previously for boridcs and other second-group elements by Schwarz ct 
al. [35] . For the H atom in hydrides the lattice contribution is even larger than 
the sphere contribution, as has been found recently [4]. It is not surprising 
that the pp component dominates the sphere contribution with the sd and pf 
components being the largest other components. 
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Table 6 

Split of the B EFGs for YB4 (optimized structure) into lattice and sphere compo- 
nents and the latter into its main contributions, i. e., sd, pp, and pf. All Vzz values 
are in units of 10^^ V/m^. 

Site Vlf V^t Vii VP/ 

B(l) -10.8 7.2 -18.0 -1.4 -14.6 -1.7 

B(2) 16.6 -4.0 20.6 1.3 17.8 1.3 

B(3) -10.8 5.0 -15.7 -1.3 -12.6 -1.6 



3.4 Knight shifts 

The Knight shifts, K, measured for both compounds, are found to be very 
small, i. e., 0.6 ppm and —1 ppm for YBg and YB4, respectively. The accu- 
racy of the K values were estimated to be ±8 ppm, where the uncertainty 
results from an incomplete reduction of the strong dipolar interaction in the 
MAS experiment, mainly between the boron nuclei. These values of K are 
much smaller than those observed in MgB2, where K(^^B) — 100 ppm [36] 
and in YNi2B2C, where K(^^B) « 400 ppm [37] (both values are for room 
temperature). This comparison indicates that the values for the densities of 
states (DOS) at the Fermi level for the boron sites are much lower than those 
in the borides mentioned above. Especially the lack of s-type DOS, which is 
usually the contribution with the largest hyperfine field, reduces significantly 
the shifts. Indeed, the s-type DOS per boron atom is, in YBg [38], over 20 
times smaller than in MgB2 [38]. Also the p-tjpe DOS per B atom is signifi- 
cantly smaller in YBg compared to MgB2. These data are compatible with the 
present electronic band structure calculations (see next Section) which show 
that s- and p-type boron states are located mainly below the Fermi energy 
and thus the contribution to the DOS at the Ep is small. 
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4 Electronic structure and chemical bonding 

The first theoretical investigation of metal borides of general formula MBg 
has been performed fifty years ago by Longuet-Higgins and Roberts [39] using 
the tight-binding approximation. According to these authors the octahedral 
arrangement of the B atoms leads to ten bonding states which, when filled 
with the 18 B s and p plus two metal valence electrons, will correspond to an 
insulator. A metal will be obtained when the hexaboride of a trivalent metal 
is considered instead of a divalent metal. In their paper the states below the 
Fermi level are not only characterized in terms of their bonding or antibond- 
ing properties but also as regards their participation in bonds within a single 
Bg octahedron or between neighbouring octahedra. Later [40] the higher ly- 
ing valence and conduction bands of YBg have been calculated by a discrete 
variational method in a non-selfconsistent fashion. Recent electronic-structure 
calculations for YBg can be found in Ref. [2] (DOS only) and in Rcf. [38] where 
also the band structure is given. For YB4 only the DOS has been published 
so far [2,38]. In both references, however, the lowest B 2s band is not shown 
although it is of crucial importance for the chemical bonding in both borides 
as will be argued below. In the case of YBg Shein et al. [38] designate it erro- 
neously as "quasi-core B 2s band" - obviously mainly because it is a low-lying 
band. In Fig. 2 the band structures for YB4 and YBq are presented. 

The DOS curves for YB4 and YBg are given in Fig. 3 together with the local 
partial DOS components. For both borides the states at energies of ca. —15 eV 
are dominated by B 2s character with a considerable amount of B 2p character. 
Such states appear for all B atoms constituting the Bg octahedra, but not for 
B(2), which belongs to the B2 unit of which each B atom is linked to two 
different octahedra. In this energy range hybride orbitals are formed which 
lead to strong covalent a bonds between the boron atoms within a Bq unit 
which can be seen from the corresponding electron densities shown in Fig. 4. 
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After an energy gap of about 3.7 eV a complex of bands appears in both 
compounds where in the lower energy range B s states prevail with an admix- 
ture of B p states except for the B(2) atoms. These states correspond mainly 
to bonds between different Eg octahedra and s-s a bonds in the B2 units of 
YB4. At higher energies below the Fermi level B p states are predominant 
apart from a relatively small B s and a considerable Y d contribution. Also 
in this energy range the interoctahedral bonds are dominant except for the 
region very near the Fermi level in YBe where the intraoctahedral B-B bonds 
again become important. 

In YBg the Y d states have mostly Cg character which can be seen from the 
valence electron densities in Fig. 5. This Figure also shows that this is not the 
case for YB4 although near the Fermi level e^-like character, i. e., {d^^, d^^-yi) 
character, is predominant. 

Our results for YBg are in very good agreement with the detailed semiem- 
pirical tight- binding investigation by Longuett-Higgins and Roberts [39]. In 
order to make the comparison we performed crystal-field splittings for the 
B p (and Y d) states. We thus obtained the distinction between their radial 
(Pz) and tangential orbitals {px,Py)- Together with electron-density plots for 
certain energy ranges, the regions for the different site orbitals for the Bg 
octahedron could be identified. We obtained the following order (increasing 
energy): aig (band at ca. —15 eV), a'lg, tiu, e'g, and t2g- (The nomenclature is 
taken from Ref. [39] and refers to the irreducible representations for the whole 
octahedron.) 

In Fig. 6 the valence electron densities in the (001) plane through the B atoms 
in YB4 and YBg is displayed. Apart from the structural differences between 
both borides (irregular heptagonal arrangement of the B atoms outside the 
octahedra in YB4 compared to an almost regular octagon in YBg) the shghtly 
higher valence electron densities between the B atoms outside the octahedra 
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in both compounds can be seen as well as the deep density minima which 
designate the much weaker covalent Y-Y bonds perpendicular to the (001) 
plane. 

5 Summciry 

For YB4 and YBg we have performed electric-field gradient and Knight shift 
measurements for the B sites. The aim of the paper has been a comparison 
of the NMR results and the available structural data with the results of accu- 
rate first-principles calculations. Wc have therefore optimized the structures 
by atomic forces and stress-tensor minimization. For the optimized structures 
EFG calculations have been performed. Perfect agreement was found between 
the experimental and calculated structural parameters and electric-field gra- 
dients thus confirming the structure models available in the literature. Based 
on calculated total and local partial DOS as well as electron-density plots, we 
analyse the bonding situation in both compounds and conclude that the low- 
est valence bands are of particular importance for the stabilty of both yttrium 
borides. 
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Fig. 1. Experimental and theoretical (model) powder-pattern ^^B spectrum for YB4 
(top) and YBe (bottom). For better readability the model spectra are shifted up 
with respect to the experimental ones. For YB4 the arrows indicate positions of 90° 
singularities (77 ~ 0) for B(2) (outermost arrows) and B(l) (middle arrows) except 
for the singularity for B(3) (innermost arrows) where the angle is different from 90° 
[rj = 0.46). 
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Fig. 2. Electronic band structure for YB4 (left) and YBg (right). 
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Fig. 3. Total DOS (full line) and local partial DOS components (B s: dotted line; 
B p: dashed line; Y d: dash-dotted line) for YB4 (top) and YBg (bottom) in units 
of states per Rydberg and per formula unit. 
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Fig. 4. Electron densities in the (001) plane through the B atoms for the lowest 
valence bands at ca. —15 eV. Left: YB4, right: YBg. A logarithmic grid of contour 
lines has been used (xj = xo2*/^). For YB4 the B(2) atoms appear as black bullets 
without a surrounding electron density. 
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Fig. 6. Valence electron densities in the (001) plane through the B atoms. Left: YB4, 
right: YBq. A logarithmic grid of contour lines has been used (xi = xq 2*/^). 
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